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We study the short time behavior of the order parameter coupled to a conserved field in semi- 
infinite geometry. The short time exponent, obtained by solving the one loop differential equations 
for the conserved density and the order parameter, agrees with the prediction from a scaling ar- 
gument based on short distance expansion. The scaling analysis further shows that this exponent 
satisfies a scaling relation similar to that known in the case of a nonconserved order parameter 
without any coupling. 
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I. INTRODUCTION 

Questions related to relaxations in critical dynamics, 
especially critical slowing down, and the nature of trans- 
port coefficients have drawn attention due to the un- 
usual properties of a system as its critical point is ap- 
proached. So far the long time relaxation has been the 
object of primary focus until recent past when the short 
time relaxation was found to reveal new universal behav- 
ior. This process sets in after the microscopic relaxation 
processes which are to be described by the microscopic 
theory. Such a short time property was first observed in 
numerical simulations and was explicitly calculated for 
a purely relaxational model with an n component non- 
conserved order parameter field <fi(x,t) (model A) ||||- 
The short time relaxation involves a new critical expo- 
nent which does not follow from any scaling relation of 
the known exponents. Being motivated by the study of 
boundary critical phenomena Q , the original study of the 
universal short time behavior was based on the consider- 
ation of a boundary in the "time" like coordinate. In the 
Renormalization Group framework there are additional 
singularities located on the time surface. In general, for 
the order parameter one may write the scaling form 0] 



-C/»o 
o 



m=(<j>(x,t))~t-^f(m t*°/S), 



(1) 



where m is the initial "magnetization" with the scal- 
ing dimension xo and £ is the dynamic exponent . 
[3 and v are the usual bulk critical exponents associ- 
ated with magnetization and length scale respectively 
j|. The function /(a) ~ a for a — > and thus for 
short time m ~ mot 8 with 8 = xo/( — P/vQ. For 
a — > oo, f(a) ~ constant and the usual long time re- 
laxation is recovered. For 2 < d < 4, since xo > f3/v, 
initially the magnetization increases upto a certain time 



t - m, 
function C(t) 
havior C(t) ~ 



In a similar fashion the autocorrelation 
= (0(x, i)</>(x, 0)) has the short time be- 
t e - d /i §. 



It is well understood now that semi-infinite systems, 
which extend over infinite space in d — 1 dimensions de- 
noted by r and over only the positive half space in the z 
direction [z > 0), have critical behavior close to the sur- 
face drastically different from the bulk Q. Detail field 
theoretical studies show that these differences arise from 
the need for an additional renormalization factor for field 
to cure the new uv singularities caused by the surface. 
Depending on the value of the surface interaction con- 
stant, conventionally denoted by c, there are different 
universality classes associated with the surface ordering. 
These are named as ordinary (c > 0), special (c = 0) and 
extraordinary (c < 0) transitions. In the ordinary transi- 
tion the surface orders along with the bulk and in the ex- 
traordinary transition the surface orders before the bulk. 
At the special point c = 0, there is a different set of expo- 
nents. The universal short time behavior is also modified 
depending on the surface universality class considered. 
As has been shown for model A ||[| , in the case of the 
special transition the order at short time grows with time 
whereas in the ordinary transition it decays with time. 
By using the short distance expansion (SDE), introduced 
by Dichl in the problem of boundary critical phenomena, 
the short time exponent in semi-infinite model A can be 
shown j^j to satisfy a scaling relation involving the bulk 
short time exponent and the static exponents. 

In this paper we are primarily concerned with the uni- 
versal short time behavior of a relaxational model with 
a conservation law. This is model C (as classified by 
Halperin et al. [Q ), where the n component order pa- 
rameter is coupled to a nonordering conserved field. We 



Ho avoid any confusion with the coordinate z, we have de- 
noted the dynamical exponent by £ instead of z 
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shall concentrate on the semi-infinite geometry. In real- 
ity such models describe binary alloys undergoing order- 
disorder transition where the concentration field (not the 
order parameter) plays the role of the conserved density 
H or systems with mobile conserved impurities or uni- 
axial Ising antiferromagnets etc. In the static limit, since 
the conserved field can be integrated out, model C be- 
comes identical to the static limit of model A with shifted 
coupling constants. The static bulk and surface exponents 
therefore can be simply borrowed from static limit of 
model A. In the dynamics of model C, the coupling of 
the order parameter with the conserved field plays an 
important role and depending on the stable fixed points 
of the parameters, the (n, d) plane can be separated into 
different regions where e.g. the relaxation rate, dynami- 
cal exponents are different 0]. The short time behavior 
for the bulk model C has been studied by Oerding and 
Janssen who, using the field theoretic renormaliza- 
tion group technique, obtain a new universal short time 
exponent (denoted as C in the following) for the order 
parameter relaxation. The dynamics of the semi-infinite 
model C has been found to be the same as the bulk 
dynamics with static exponents same as model A at dif- 
ferent universality classes of surface transitions. 

The short time dynamics and its universal features are 
important in quenching experiments where the system is 
taken from an unstable ordered state at high tempera- 
ture to, say, the critical temperature. In numerical sim- 
ulations ]l2] ] , the short time dynamics is relatively easier 
to observe because in this time regime the critical slowing 
down does not set in. In several cases the results from 
early time dynamics gave good estimates for not only 
the short time but also the bulk exponents. Therefore 
the knowledge about the short time exponent for model 
C and its relation to other critical exponents is expected 
to be useful for numerical and experimental work on anti- 
ferromagnetic systems [O, binary alloys and other 
systems where there are coupling between the order pa- 
rameter and a conserved density field. 

A priori it is not clear what role a conserved density 
would play in the short time dynamics of the semi-infinite 
model C. Another related question is whether the short 
time dynamics can be explained by SDE in a similar 
manner as model A. By solving one loop equations of 
motion for the conserved density and the order param- 
eter, we show that there are certain subtleties in this 
situation. Close to the surface the conserved density has 
spatial variation not coming from SDE. However, this 
contributes to the order parameter equation significantly 
and finally leads to the order parameter relaxation con- 
sistent with SDE. Another interesting goal of studying 
model C is to observe the shape of the conserved den- 
sity profile. This question is partially answered in our 
analysis in a region very close to the surface. 

The dynamics of the order parameter field and the con- 
served nonordering field follows from the hamiltonian 



H[<j>, E] = J d d x{ 



^ 2 +i^+^+i E <p 2 ) ( 2 ) 



and the Langevin equations 



<9 t 0(x,i) = — A- 



SH 



(x,t) 



d t E(x,t) = XpV 2 ( 



sn 



SE(x, t) 



+ C(*,t) 



(3) 
(4) 



The Gaussian random noise 77(2, t) and C,{x,i) have zero 
mean and correlations 

(C(x, f)C(x', t')) = 2A S(x - x')6(t - t') (5) 
(r/(x, i)?7(x', t 1 )) = -2A p V 2 <5(x - x')5(t - t'). (6) 

The generating functional in terms of the response 
fields 4>(x,t) and £(x,t) is (TsUll 



J[<j), j>, m, fh] = dt d d x[</>{d t( />(x, t) + A(r - V 2 )</> + 



^ 3 + A 7 <Kx,i)£(x,<)} - A0 2 (x, t)+E(x, t)d t E{x,i) 



\p(V 2 E(x, t))( j0(x, tf + E( X , t)) - \p(VE) 2 } 



(7) 



The correlation propagator in the bulk case contains an 
equilibrium part which is translationally invariant in time 
and a non equilibrium mirror symmetric part. To be 
brief, we mention only the correlation propagator for the 
field <j> and refer the reader to |10| for the other response 
and correlation propagators. In Fourier space the corre- 
lation propagator is 

CV(q,i-0 = <<Kq,i)<H-q,0) = [e^^l^'l - 

e -Kr+ q 2 ){t+t')y {T + q 2 ) (8) 

This propagator corresponds to the Dirichlct initial con- 
dition. 

The static version is equivalent to the static limit of 
model A with <p A coupling u = g — 3-f 2 . We briefly recall, 
from the renormalization group analysis of the dynamics 
of model C Q, the results useful for the present work. 
The dimensionless coupling constants u and v are de- 
fined as KdU = p e u and K^j 2 = p e v, where p is an 
arbitrary momentum scale and e = 4 — d and K^ 1 — 
2 d ~ 1 7r d / 2 r[d/2]. Broadly there are two distinct regimes 
in the (n, d) plane depending upon the sign of the specific 
heat exponent a. For a < the stable fixed points are 
u* = 6e/(n + 8) and v* — 0. Clearly in this regime large 
scale properties of model C are same as model A and the 
dynamic exponent ( = 2 + 0(e 2 ). For a > there is a 
stable nonzero fixed point v* = 2e(4 — n)/[n(n+8)]. These 
two regimes can further be separated depending on the 
fixed point of p. For the details on these aspects we refer 
the reader to For one component order parameter 
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p* = 1 + 0(e) and the conventional scaling £ = 2 + ajv 
holds good. 

In an attempt to understand the short time exponent 
for model C in a semi-infinite geometry using a scaling 
argument based on a SDE, we write the scaling form for 
the magnetization 



m(z,t,m ) ~ r^ffz/t'/^mot 10 ^). 



(9) 



Since the conserved density profile has a nontrivial be- 
havior solely due to the coupling with the order parame- 
ter field, we have ignored the explicit dependence of the 
scaling function on the conserved field. Our approach of 
solving one loop renormalized equation of motion pro- 
vides a more rigorous justification for this. Furthermore 
because of the scaling of the magnetization close to the 
surface asm~ t-Pi/"( ) where /3i is the surface magne- 
tization exponent |J], we expect !F(x,y) ~ x a J-\(y) such 
that a = —f3)/v. This implies that near the boundary 



(10) 



This result can also be understood by SDE which should 
hold good for fields with different scaling dimensions on 
or off the surface. Since at short time m should be pro- 
portional to mo, the short time behavior close to the 
surface is 



m 



(11) 



Eqn. ( |ll| ) shows that the early time behavior close to the 
surface is described by the exponent B c \ which satisfies 
the scaling relation 



09 - /3i)M, 



(12) 



involving the bulk short time exponent and static expo- 
nents. 

The above scaling relation is supported by an explicit 
calculation starting from the linearized equations for the 
order parameter and the conserved density. This ap- 
proach is elaborated in section II. Singularities appearing 
at the one loop level are taken care of by the renormaliza- 
tion of appropriate parameters. Solving the renormalized 
equations, the short time exponents are then obtained. 
The details in Section II is important to appreciate the 
crucial role played by the conserved density. The justi- 
fication of the scaling argument and the agreement be- 
tween the two approaches are discussed in section III. 



II. ONE LOOP RENORMALIZED EQUATIONS 

Starting from initial conditions and the translational 
invariance in d — 1 directions, the one loop equations 
for time variations of the averaged conserved density 
£{z,t) = (E(x, t)) and order parameter m(z,t) = 
(</>(x, t)) can be written as 



1 

V 
1 

r 



— n l S = d 2 z £(z,t) + ^d 2 z (<p 2 (x,t)) 

■ i / i \ I r ^2 , u {n + 2) 
—OtTn{z, t) + [Tom — V m H 



(13) 



C{z,t)m + 



l0 (E(z,t)cb(z,t))}=0, 



(14) 



The one loop term introduces spatial variation in the 
conserved density profile which is flat otherwise. The one 
loop contribution C(z, t) — (</> 2 (x, i)) consists of bulk and 
surface parts as C(z,t) = C&(0,i) ± Cb{2z,t) |17|, where 



C b (z,t) 



(2ir) d J r a + q 



M _ e -2A(g 2 +r)t] 



exp[ifcz] 



(15) 



and ± refer to special and ordinary transitions respec- 
tively. The first term on the RHS of (|l5|) in the bulk 
contribution Cb(0, t) has a divergence which has to be ab- 
sorbed by the renormalization of the temperature. Eval- 
uating the integrals we obtain 



C(z,t) 



32tt 



' [4±^exp(-z 2 /2t)], 



(16) 



where t = Xt. Substituting this in (|l3|), we find that the 
solution of the conserved density profile is of the form 
£{z,t) — =F(z 2 /t), where t = Xpt. Very close to the 
boundary such that z <C t 1 ^ 2 , F(x) ~ T 32^30 where — 
refers to the special transition. Therefore very close to 
the boundary we have 



£( z ,t) - T7/32tt 2 



(17) 



for special and ordinary transition respectively. Eqn. 
(|l7| ) shows the importance of the one loop term from <p 
(the inhomogeneous term in (|l3[)) in the behavior of (E) 
near the surface. Since the above result is restricted to 
the regime z i 1 / 2 , initial condition cannot be reached 
from this. At the fixed point of our interest this power 
law form associated with a prefactor (^(e 1 / 2 ) contributes 
at 0(e) in the equation for the order parameter. Though 
a further analysis about the shape of the conserved den- 
sity profile away from the surface deserves attention Q , 
we here restrict ourselves very close to the surface. 

Next we consider the term A(x,t) = (E(x,t)4>(x,t)) 
which needs to be expanded in order to take into ac- 



count the other O(e) terms in (14). Using the generating 
functional in (Q), we have 

A{x,t) = (E{x,t)(f)(-xL,t)) -A7 J d d Xl J dti 

{(^(x,<) J E(x 1 ,t 1 )}(0(x,i)0(x 1 ,t 1 )}(0(x 1 ,t 1 ))- 
p<S(x,i)V 2 ^(x 1 ,t 1 ))( ( /»(x,i)0(x 1 ,t 1 )}(0(x 1 ,t 1 )}} (18) 

For convenience we denote the two terms in the curly 
bracket by A\(x, t) and A%(x, t) respectively. It is appar- 
ent that there is no straight forward way to evaluate the 
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last two terms due to their coupled structure. To obtain 
the contributions of these term in the bulk case, we as- 
sume (0(xi,ti)) to be space time independent. This is 
also justified if, say, the time dependence has an univer- 
sal exponent ~ t p where p ~ 0(e). In that case from 
the expansion (0(xi,ti)) ~ 1 + e\nt it is clear that upto 
this order of calculation only the time independent piece 
is the important one. Now the contributions from the 
above terms in the bulk case are 



Ax(x,t) 



s d 7 , 



1 



1 



A 2 (x,t) = - 



(2ir) d p2(l + p) 2 4p(l + p) J 
Sd 7, P , Pi 



(27r) d p4(l + p) 2(1 + p) 2 



(19) 



where Sd — r[d'/2 ] • These two terms correspond to second 
and third diagrams in Fig 1 of jfij. Correspondingly the 
second terms in the square brackets can be compared 
with 0(1) parts in the fourth and third terms in Eqn. 
(32) of pel . These are the contributions from the 'ini- 
tial' part of the correlators. 

In the semi-infinite case we have 



A 1 (z,t) = ± 



s rf 7 2 j oexp[-z 2 /(l + p)t] exp[—z 2 /2pi\ 



A 2 (z,t) = ± 



(27r) dL 2z 2 {p 2 -\) 2z 2 (l-p) ' 

s d p7 2exp[-z 2 /(l + p)t] exp[-z 2 /2i]. 



(2ir)° 



2z 2 {1 - p 2 ) 



2z 2 {\- p) 



Adding all the bulk and surface contributions from above 
we have 



A 



dtm — d 2 m 



(64tt 2 



-l± 
t 



\m 



7 



IGttH 
37 2 to 



-TO 



=F 



32tt 2 z 2 



= 



(20) 



Note that for 7 = 0, we get back the equation for 
model A. This equation can be solved by assuming a 
scaling form m (z,t) = UfyViz/t 1 ^). Here we are re- 
stricted to £ = 2. In the bulk limit z — ► 00 we obtain 
U(t) ~ P", where 9 C = 0{( 2 ). This agrees with the 
bulk short time exponent in |u|. For z <C i 1 / 2 , we find 
V(z/t 1 ^ 2 ) ~ (z/i 1 / 2 ) 1-6 / 6 for the ordinary transition and 
Viz/t 1 / 2 ) ~ (z/t 1 / 2 )-^ 6 for the special transition. Thus 
for magnetization upto 0(e), we have 

m(z, i)~ m z 1 ~ e/e t~ 1/2+e/12 for ordinary transition (21) 
~ TOoz _e ^ 6 i £ ^ 12 for special transition. (22) 



III. CONFORMITY WITH SCALING ANALYSIS 

A few points in the above perturbative calculation and 
in the previous scaling analysis need to be re-emphasized 
here. From the nature of the one loop equation for the 
conserved density and its solution, it is clear that we are 



in a time regime where the conserved density profile is 
controlled by the <f> correlations. As a consequence, the 
inhomogeneous equation ( |l3| ) provides a spatially depen- 
dent one loop correction to the conserved density profile 
E with a prefactor of 0(e 1 / 2 ). The conserved field is, 
therefore, redundant in the scaling function which in- 
volves only scaled variables with e dependence in various 
exponents. 

We see that the short time exponents in (^l|) and 
(p2|), obtained by solving the differential equations, are in 
agreement with the prediction from scaling analysis upto 
the factor £ (recall that for bulk model C, 9 C = 0(e 2 ) 
p0{ and the static exponents (3 and (3\ are same as ob- 
tained from the static limit of model A). Clearly while 
solving one loop equations, we are restricted to ( — 2 
because of the very nature of the equation. The fact 
that C = 2 + e/3 should appear from the appropri- 
ate propagator renormalization which has not been per- 
formed in this simple approach. Since the static expo- 
nent of z, obtained from the scaling analysis, is in clear 
agreement with the calculation, we use ( — 2 + e/3 in 
V{z/t 1 ^) - (z/t 1 ^)-^ 6 , (2r/t x /C)i-e/6 for S p ec i a i and 
ordinary transitions respectively. Thus we conclude that 
the short time exponents in the semi-infinite model C are 

6 c i= e/6C for special transition (23) 
= —(1 — e/6)/C for ordinary transition. (24) 

The characteristics of model C is reflected in the expo- 
nent of time, whereas the exponent of z remains same as 
model A. This is due to the fact that the static limits of 
model C and model A are same. We do not discuss more 
about the behavior of the autocorrelation function. The 
short distance, short time behavior of this quantity for 
model A has been explained in B. Our analysis shows 
that the same scaling relation is valid here. 



IV. CONCLUSION 

We have obtained the nature of the conserved density 
profile close to the surface by solving a one loop differen- 
tial equation. This one loop correction to the flat profile 
is of 0(e 1 / 2 ). This form of the profile has significant 
impact on the relaxation of the order parameter in the 
semi-infinite geometry. The universal short time expo- 
nent 9 c i in a semi-infinite geometry is given in (p3[ ) and 
(pi|). As in model A the short time exponents for spe- 
cial and ordinary transitions differ drastically and they 
are consistent with the scaling analysis based on short 
distance expansion. 
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